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JE] 要 : 针对 维度 学 习 策 略 (Dimensional learning strategy, DURP qe) t 度 开发 ”问题 ， 提 出 了 一 种 综合 维度 学 习 
的 多 群 协作 粒子 群 优化 算法 (CDL-MCPSO)， 为 提高 种 群 搜索 效率 ， uU MP A 雪 构 ， 将 种 群 划分 
为 一 个 主 群 和 四 个 从 群 ， 主 群 执行 综合 学 习 策 略 在 搜索 空间 进行 大 范围 探索 ， 从 群 执行 综合 维度 学 习 策 略 
(Comprehensive Dimensional Learning，CDL) 在 局 ee 度 地 开发 ， 主 从 群 通过 执行 具有 不 同 职能 能 的 
算法 能 够 有 效 实现 其 在 勘探 和 开发 之 间 的 平衡 ， 同 时 为 保持 种 群 多 样 性 ， 提 出 了 一 种 新 的 解 交 换 机 制 (SEM)， 用 来 
在 主 从 群 独立 运行 各 自 算 法 若干 代 之 后 进行 信息 的 交流 与 协作 ， 以 指导 粒子 后 期 进行 更 准确 的 搜索 ， 最 后 ， 针 对 初 
始 化 过 程 随机 性 过 高 ， 运 用 拉丁 超 立 方 体 采样 方法 对 算法 重建 输入 分 布 ， 为 验证 CDL-MCPSO 的 有 效 性 ， 将 其 与 5 
种 粒子 群 算法 变 体 在 10 个 测试 函数 进行 实验 对 比 ， 结 果 表 明 该 算法 总 是 可 以 找到 优 于 或 相当 于 对 比 算 法 的 解 ， 在 
求解 复杂 函数 时 具有 可 行 性 和 高 效 性 
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Multi-swarm collaborative particle swarm optimization algorithm 
based on comprehensive dimensional learning 


Zhang Qiwen, Wang Yangting! 
(School of Computer & Communication, Lanzhou University of Technology, Lanzhou 730050, China) 


Abstract: Aiming at the problem of "over-exploitation" in dimensional learning strategy (DLS) , this paper proposed a multi- 
swarm collaborative particle swarm optimization algorithm with comprehensive dimension learning (CDL-MCPSO) . The 


cluster structure of master-slave paradigm divides the population into a master group and four slave groups, in which the 
master group executes a comprehensive learning strategy to conduct a large-scale exploration in the search space, and the 
slave groups execute a comprehensive dimensional learning strategy (CDL) to exploit high-precision solutions near the local 
optimal. The master-slave groups can effectively achieve the balance between exploration and exploitation by executing 


algorithms with different functions. At the same time, in order to maintain the diversity of the population, the CDL-MCPSO 
proposed a new solution exchange mechanism(SEM) , By doing this, the algorithm can accomplish information exchange and 
cooperation after the master-slave groups running their respective algorithms for several generations independently, so as to 
guide the particles to conduct more accurate searches in the later stage. Finally, for the high randomness of the initialization 
process, the algorithm adopted latin hypercube sampling to reconstruct the input distribution. In order to verify the 
effectiveness of CDL-MCPSO, compared with 5 kinds of PSO variants in 10 test functions. the results show that the algorithm 
can always find better or equivalent solutions. It is feasible and efficient in solving complex functions. 

Key words: particle swarm algorithm(PSO); comprehensive dimensional learning strategy; master-slave paradigm; multi- 
group collaborative 


0 ”引言 算法 性 能 的 提升 有 很 大 作用 ， 进 而 在 基本 粒子 群 算法 中 引入 

一 了 惯性 权重 ，Sarhani 等 人 BJ 为 平衡 算法 的 全 局 与 局 部 搜索 能 
近 些 年 来 ， 智 能 优化 算法 逐渐 成 为 解决 科研 和 工程 应 用 力 ， 提 出 一 种 基于 隐 性 马尔 可 夫 模 型 状态 估计 的 惯性 权重 控 
领域 中 许多 复杂 优化 问题 的 首选 方法 。 常 见 的 主要 包括 粒子 。 制 机 制 ， 实 验 结果 证 明了 算法 的 收敛 优势 。2017 年 ，Gou 等 
群 算法 (PSO)、 和 声 搜索 算法 (HS)、 文 化 算法 (CA)、 基 于 教 与 人 外 根据 每 个 粒子 之 间 的 个 体 差异 动态 调整 参数 的 值 ， 同 时 
学 的 算法 (TLBO) 以 及 重力 搜索 算法 (GSA) 等 。 其 中 最 经 典 的 。” 采用 重启 策略 再 生 粒 子 ， 实 验 结果 证 明了 该 算法 具有 和 鲁 棒 性 
当 属 模拟 生物 行为 机 制 的 粒子 群 优化 算法 (particle warm ”和 可 扩展 性 。 


ih 


optimization, PSO)!!, 'Ziégi Kennedy 和 Eberhart 博士 于 在 拓扑 结构 的 选择 方面 , Nasir 等 人 J 提出 了 一 种 动态 邻 

1995 年 正式 提出 的 一 种 基于 种 群 的 元 启发 式 优化 算法 ,由 于 域 学 习 粒 子 群 优 化 算法 (DNLPSO)， 该 算法 从 邻 域 (包括 自身 ) 

其 实现 简单 、 控 制 参数 少 吸 引 了 大 量 学 者 研究 并 涌现 出 了 许 的 最 佳 位 置 选 择 一 个 样本 粒子 。 文 献 [6] 为 了 从 根源 上 克服 

多 变 体 。 针 对 它 的 改进 主要 有 四 个 方面 : 参数 设置 、 拓 扑 结 PSO 算法 易 陷 入 局 部 最 优 的 问题 ， 针 对 多 种 拓扑 结构 进行 实 

构 、 算 法 融合 以 及 学 习 策 略 。 验 分 析 ， 最 终 提出 了 一 种 全 局 探索 与 局 部 开发 相 融 合 的 混合 
参数 的 取 值 对 于 平衡 算法 的 勘探 和 开发 能 力 有 着 重要 的 拓扑 粒子 群 算法 。 

EH o Shi 等 人 咏 通 过 研究 发 现 粒 子 对 先前 速度 的 继承 能 力 对 利用 不 同 的 进化 算法 来 提高 PSO 的 性 能 是 研究 人 员 的 
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重点 。 最 早 的 算法 融合 是 Gong 等 人 忆 将 遗传 算法 中 的 “能 力 变 成 之 前 的 两 倍 ， 种 群 多 样 性 逐渐 丧失 。 针 对 此 不 足 ， 
三 个 经 典 操作 选择 、 变 异 和 交叉 引入 到 PSO 中 。 每 个 粒子 根 CDL-MCPSO 通过 融合 综合 学 习 策 略 (CLS) 和 维度 学 习 策 略 
据 遗 传 操作 生成 一 个 有 希望 的 样本 , 然后 以 标准 PSO 的 方式 (DLS) 各 自 的 优点 ， 提 出 “综合 维度 学 习 策略 (CDL)”, 去 掉 公 
从 样本 中 学 习 。2018 Æ, Aydilek 等 人 外 将 萤火虫 算法 与 粒子 。” 式 第 三 部 分 对 全 局 最 优 项 的 学 习 , 引入 CLS 中 对 所 有 其 他 粒 
群 算法 的 概念 进行 融合 ， 实 验 结果 证 明了 算法 不 仅 能 有 效 平 ” 子 概率 性 学 习 的 思想 ,使 粒子 在 开发 的 同时 也 可 以 适当 探索 ， 
衡 算法 的 探索 与 开发 能 力 ， 而 且 运行 时 间 和 解 的 精度 都 有 明 ”并 基于 信息 流动 原理 提出 了 一 种 新 的 解 交 换 机 制 (SEM)， 在 
显 提升 。Senel 等 人 中 将 PSO 的 开发 能 力 与 灰 狼 优化 算法 的 是 高 算法 收敛 速度 的 同时 也 优化 了 解 的 质量 ,最 后 通过 拉丁 
探索 能 力 相 结 合 ， 在 同一 测试 环境 下 ， 该 算法 在 短 时 间 内 获 。 超 立 方 体 采样 对 算法 重建 输入 分 布 ， 有 效 避 免 了 初始 化 过 程 


得 较 好 的 成 效 。 实 验 结果 表明 算法 能 在 更 短 的 时 间 内 找到 更 于 随机 性 过 高 给 算法 带 来 的 负面 影响 。 通 过 对 改进 思路 / 
好 的 优化 方案 。 生 原 因 的 描述 ,可 以 看 出 CDL-MCPSO 具备 新 颖 性 和 前 沿 性 。 


根据 以 上 各 种 改进 , 可 以 明确 PSO 变 体 的 主要 目标 是 保 1 ”相关 工作 
持 种 群 勘探 和 开发 能 力 的 平衡 ， 并 在 保持 种 群 多样 性 的 同时 
提高 算法 的 收敛 速度 以 及 解 的 质量 ， 但 现实 生活 中 大 多 数 优 ”1.1 粒子 群 优化 算法 
化 问题 往往 都 比较 复杂 。 在 多 样 性 方面 ， 单 维 问题 需要 较 小 PSOP 源 于 对 鸟 群 员 食 共享 个 体 知识 社会 行为 的 研究 。 
的 种 群 多 样 性 ， 以 便 快速 收 僵 ， 而 多 模 态 问题 需要 较 好 的 种 。 群 中 的 每 个 个 体 称 为 粒子 ， xs 个 潜在 解 ， 
群 多 样 性 来 避免 算法 陷入 局 部 极 值 。 由 此 可 见 只 采用 单一 策 ”粒子 通过 跟随 自身 个 体 最 优 和 种 群 全 局 最 优 来 更 新 自己 的 速 
略 无 法 满足 不 同 优化 问题 的 需求 。 因 此 一 些 研究 人 员 提 出 了 度 和 位 置 ， 直 到 得 到 满足 终止 条 件 的 最 优 解 。 粒 子 速度 和 位 


分 而 治之 的 策略 来 解决 复杂 优化 问题 ， 针 对 不 同 的 优化 目标 置 更 新 公式 如 下 : 

运用 的 不 同 的 学 习 策略 。 陈 跃 刚 等 人 09 基 于 综合 学 习 策略 对 Vj(t+D= VD +on (CI (0) - X; O+ " 
多 目标 问题 进行 求解 ,通过 分 解 方法 更 新 主导 粒子 以 增强 解 con, QC" CD 一 XGOD) 

的 分 布 ;采用 存档 机 制 存 储 优化 过 程 中 的 非 支 配 解 ,并 采用 多 X, (t€) 5 V. (t+D+ X40 (2) 


项 式 变 异 来 避免 陷入 局 部 最 优 , Mendes 等 人 00 提 出 了 完全 知 其 中 :i=142.…m,j=42..…n, m 和 分别 是 种 群 规模 和 维 数 。@ 
情 粒 子 群 优化 算法 (FIPS)。FIPS 利用 完全 联通 邻 域 中 所 有 邻 RERE, a e 分 别 是 自我 认 知 和 社会 i 
居 的 个 人 最 佳 信息 的 平均 值 更 新 粒子 的 位 置 ， 文 献 [12] 在 考 ”第 了 维 上 [0,1] 之 间 均 匀 分 布 的 两 个 随机 数 。 XE RI XI 分 别 
虑 了 种 群 结构 、 多 模式 学 习 和 个 体 间 博 弈 等 因素 之 后 ， 提 出 民 表 粒子 个 体 最 优 和 全 局 最 优 。 
了 具有 博弈 概率 选择 的 多 子 群 粒子 群 算法 。Parsopoulos 等 人 器 1.2 综合 学 习 策略 (CLS) 
是 出 了 统一 粒子 群 优化 算法 (UPSO), 通过 选择 合理 的 邻 域 大 在 2006 年 Liang 等 人 P0 提 出 的 CLPSO 采用 新 的 学 习 策 
小 ， 利 用 全 局 邻 域 和 局 部 邻 域 的 最 佳 信息 来 更 新 粒子 位 置 ， 各 ， 通 过 利用 所 有 其 他 粒子 的 个 体 最 优 信 息 来 更 新 速度 ， 为 
囊 小 平等 人 (4 将 种 群 动态 地 划分 为 三 个 不 同 阶 层 ， 并 根据 不 粒子 提供 了 更 多 的 学 习 标本 和 潜在 的 搜索 空间 ， 在 综合 学 习 
同 阶层 粒子 特性 ,分 别 采用 局 部 学 习 模 型 、 标 准 学 习 模 型 以 及 ”策略 中 , 粒子 的 更 新 由 学 习 样本 X% 四 来 指导 ,速度 更 新 公式 


全 局 学 习 模 型 , 来 改善 算法 性 能 , 文献 [15] 开 发 了 两 种 新 的 学 。 ”如 下 : 

习 策略 ， 一 方面 ， 提 出 了 适应 度 景观 中 的 局 部 稀 疏 度 度量 方 VD - OV D+er XsD - X.) Q) 
法 来 估计 粒子 的 拥塞 和 分 布 情况 ， 并 在 此 基础 上 通过 引导 粒 位 置 更 新 见 式 (2)， 其 中 X50 表示 粒子 ;个 体 最 优 pbest 
子 到 稀疏 区 域 来 建立 探测 策略 。 另 一 方面 ， 基 于 多 群 策 略 和 ” 的 第 J 维 通过 向 所 有 其 他 粒子 个 体 最 优 的 第 7 维 进行 学 习 
自 适应 子 群 大 小 调整 提出 了 一 种 自 适应 开发 策略 。 Zhou 等 人 09 (包括 粒子 i 自身 个 体 最 优 的 第 j 维 ) 后 产生 的 样本 ， 样 本 
提出 了 一 种 自 适 应 分 层 更 新 粒子 群 优化 算法 ， 为 两 个 群 分 别 X50 学 习 过 程 : 粒子 ;个 体 最 优 每 一 维 具体 向 谁 学 习 取决 于 
生成 两 层 和 三 层 更 新 公式 ,并 引入 了 多 项 选择 综合 学 习 策略 。 ”学 习 概 率 p.， 不 同 的 粒子 学 习 概率 不 同 ， 首 先 为 每 一 维 生 成 


Şaban 等 人 "1 基于 综合 学 习 和 主 从 范式 提出 了 并 行 综 合 学 习 ” ”一 个 随机 数 rand， 若 随机 数 大 于 p. ， 则 向 自己 个 体 最 优 相应 
粒子 群 优化 算法 (PCLPSO)， 它 基于 主 从 范式 ， 并 在 不 同 范式 ”维度 学 习 ， 否 则 向 其 他 粒子 个 体 最 优 相 应 维度 学 习 ， 学 习 概 


> 间 进 行 解 信息 的 协作 。 文 献 [18] 受 天 体 学 和 植物 学 启发 , 提 ”” 率 p. 更 新 公式 如 下 : 

出 一 种 多 策略 融合 的 粒子 群 优化 算法 。 其 中 三 黑洞 系统 捕获 (exp 7D, 

策略 和 多 维 随机 干扰 策略 用 来 增强 粒子 的 全 局 搜索 能 力 ， 而 P. =0.05+0.45x ps-l (4) 

协调 因子 用 来 完成 从 全 局 寻 优 向 局 部 搜索 的 转变 ,以 提高 收敛 速度 。 | exp(10) -1 

文献 [19] 根 据 PSO 算法 的 收敛 特性 和 Logistic 映射 的 混沌 思想 , 提 ps 代表 种 群 规模 ， 在 综合 学 习 策 略 中 ， 其 他 粒子 先前 的 

出 一 种 粒子 置换 的 双 种 群 综合 学 习 PSO 算法 (PP-CLPSO)。 最 佳 位 置 是 任何 粒子 可 以 学 习 的 样本 ， 并 且 粒 子 的 每 个 维度 
从 以 上 改进 方法 不 难看 出 ,大 多 数 针 对 参数 、 拓 扑 结构 、 都 可 以 从 不 同 的 样本 中 学 习 ， 这 种 融合 来 自 不 同学 习 样 本 信 

算法 融合 的 改进 ， 仅 仅 是 针对 其 中 一 两 个 目标 进行 优化 , 没 ” 息 的 策略 扩大 了 粒子 的 潜在 搜索 空间 ， 大 大 提高 了 种 群 的 全 


有 平衡 好 多 个 目标 之 间 的 关系 ， 多 种 策略 分 而 治之 思想 是 根 局 搜索 能 力 ， 因 此 在 解决 多 模 态 问题 效果 优异 ， 但 由 于 去 掉 
据 算 法 的 不 同 的 优化 目标 设计 不 同 的 学 习 策略 ， 不 仅 能 够 实 。 了 对 全 局 最 优 项 的 学 习 ， 导 致 在 求解 单 峰 问 题 上 效果 不 佳 。 
现 勘 探 和 开发 的 平衡 ， 还 能 有 效 避 免 早熟 收敛 、 提 高 收敛 速 1.3 维度 学 习 策略 (DLS) 


度 以 及 解 的 质量 。 Xu 等 人 RWI 针对 CLPSO 不 足以 及 标准 PSO 在 种 群 更 新 
CDL-MCPSO 和 以 上 多 策略 融合 的 改进 一 致 ， 也 采用 中 存在 的 “振荡 ”P22 和 “两 步 走 ， 一 步 回 头 ”P3 问 题 ， 提 出 


分 而 治之 的 思想 ,但 是 考虑 了 绝 大 多 数 算法 没有 考虑 到 的 “ 振 “了 维度 学 习 策 略 ， 利 用 构造 的 学 习 样 本 XI 和 全 局 最 优 x em 
荡 ” 和 “维度 退化 ”现象 , 例如 文献 [20] 提 出 的 维度 学 习 策略 ， 更 新 速度 , 学 习 样 本 x? 的 构造 : 区 别 于 综合 学 习 策略 中 学 习 
虽然 解决 了 “维度 退化 ”现象 ， a a 样本 Xe 的 每 一 维 来 自 其 自身 或 所 有 其 他 粒子 的 相应 维度 

学 习 策 略 并 不 完全 适 配 子 群 职能 而 产生 了 新 的 问题 -“ 开 采 过 DLS 学 习 样 本 x? 的 每 一 维 是 来 自 其 自身 或 全 局 最 优 
E", 原因 在 于 公式 第 二 项 的 学 习 样 本 是 全 局 最 优 , 在 少数 迭 xee 的 相应 维度 ， 学 习 与 否 取决 于 维度 替换 后 粒子 适应 度 值 
代 之 后 ， 会 和 公式 第 三 项 完全 一 致 共同 发 挥 开 发 作用 ， 开 发 ”是 否 变 好 。 这 种 策略 可 以 将 x o 的 优秀 信息 传递 给 样本 XP, 
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录用 定稿 
促进 了 
公式 如 下 : 


Xm 编码 模式 的 传播 和 保护 ，。 


张 其 文 ， 等 : 


Vij(t+D= oV, Oan (X5 0)- X,,0)* 


Ch, 


APTO- X0) 


位 置 更 新 按 
体 最 优 xr” 整合 
学 习 后 的 样本 不 
学 习 项 添加 了 回 
也 大 大 加 快 了 算 


代 之 后 相差 不 大 ， 


照 标 准 粒 子 群 
全 局 最 优 x ee 
会 比 个 体 最 亿 
来 , 不 仅 解决 了 
法 的 收敛 速度 ， 


有 前 途 维 
ti 差 ， 同 时 将 对 全 局 


更 新 式 (2), 维度 学 习 策略 让 个 


维度 学 习 策略 速度 更 新 


G) 


度 上 的 信息 


， 保 证 了 
fg fü xe 的 


“振荡 ”和 “维度 退化 ”现象 ， 


但 由 于 公式 后 


项 在 更 新 几 
均 侧 重 于 开发 ， 极 易 出 现 “ 过 度 开 发 


» 
o 


2 综合 维度 的 多 群 协作 粒子 群 优化 算法 (CDL-MCPSO) 


UEZ E 


习 策 略 


位 置 在 更 新 过 程 


于 忽略 了 维 


中 会 出 现 “前 进 两 步 ， 


即使 粒子 适应 度 


全 局 最 优 的 程度 大 于 其 他 维度 
导致 粒子 在 靠近 和 远离 全 


度 退 化 现象 ， 


使 得 粒子 


于 退 一 步 ” 的 问题 ， 


度 变 好 了 ， 


8 A xD 为 在 3 些 维度 


度 上 粒子 靠近 
上 粒子 远离 全 局 最 优 的 程度 ， 


综合 维度 学 习 的 多 群 协作 粒子 群 优 化 算法 
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子 个 体 最 优 对 其 他 所 有 粒子 个 体 最 优 的 学 习 ， 使 得 到 的 CDL 
可 以 适当 削弱 DLS 的 开发 程度 ， 避 免 粒子 过 早 陷入 局 部 最 
优 ， 发 挥 出 DLS 最 好 的 优化 效果 。 速 度 更 新 公式 具体 如 下 : 

VC+D=aV (0) +on OC) 7 X0) + 

on (XL) - X; (0) (6) 

位 置 更 新 见 式 (2)， 所 提 4 策略 是 对 标准 粒子 群 算法 更 新 公 

式 第 二 项 和 第 三 项 换 种 思维 的 改进 ， 可 以 使 算法 在 具有 
开发 能 力 的 前 提 下 适度 探索 ， 由 于 CDL 学 习 样 本 对 粒子 
索 方向 的 强大 指引 ， 整 体 上 开发 能 力 更 胜 一 筹 ， 因 此 在 集群 
结构 中 用 来 对 从 群 进行 更 新 ， 主 群 负责 在 大 范围 搜索 空间 进 
行 探索 ， 因 此 选择 全 局 搜索 能 力 强 的 CLS， 具 体 更 新 过 程 见 
1.2 小 节 。 
2.3 ”多 群 协作 
Vac cu eU E 
量 的 解 ， 而 且 还 能 很 好 地 保持 种 群 多 样 性 ， 会 使 算法 的 
全 速度 变 慢 ， 从 群 负责 在 较 好 解 附近 进行 高 精度 的 开发 


E 


局 最 优 之 间 反 复 横 跳 ， 降 低 了 算法 能够 大 大 加 快 算 法 收敛 速度 ， 但 是 若 此 时 靠近 的 全 局 最 优 是 
搜索 效率 也 影响 了 解 的 质量 , 同时 本 文 在 对 DLS 速度 更 新 公 ”局 部 最 优 ， 会 使 粒子 无 法 逃离 此 区 域 ， 综 上 分 析 ， 需 要 利 
式 演算 之 后 发 现 , 公式 第 二 项 的 XUS 在 向 Xs" 学习 几 代 后 构 。” 用 从 群 搜索 到 的 信息 来 指导 主 群 加 快 收敛 速度 ， 也 需 利 用 
造 的 学 习 样 本 x? 会 慢 慢 地 非常 接近 全 局 最 优 Xw* 的 值 ， 此 ” 主 群 搜索 的 信息 来 帮助 从 群 粒子 逃离 “ 局 部 最 优 ”。 故 主 从 
时 对 于 后 面 粒子 的 更 新 ， 公 式 后 两 项 的 学 习 样本 均 为 全 局 最 。 群 之 间 进 行 信息 交流 是 必要 的 ， 交流 过 程 需要 考虑 以 下 四 
优 ， 赋 予 公式 最 后 一 项 的 开发 能 力 由 于 后 两 项 的 无 限 接近 变 。 个 问题 :(1) 迁 移 周 期 ，(2) 交 流 的 对 象 ，(3) 需 要 交换 的 信息 ;(4) 整 
成 了 两 倍 ， 而 且 粒子 在 更 新 过 程 中 没有 向 其 他 粒子 学 习 ， 没 合 策略 。 为 解决 以 上 问题 ， 提 出 了 一 种 新 的 解 交 换 机 制 (SEMD。 
有 可 以 帮助 粒子 逃离 “局 部 最 优 ” 的 条 件 ， 极 易 出 现 “ 开 发 (1) 迁 移 周 期 。 表征 的 是 主 从 群 在 独立 执行 各 自 算法 多 少 
过 度 ” 为 解决 此 问题 并 提升 算法 运行 效率 , 本 文 首先 采用 主 。 代 之 后 开始 进行 信息 的 交流 与 协作 ， 不 管 什 么 算法 在 优化 一 
从 范式 的 集群 结构 :一 个 主 群 和 四 个 从 群 ， 其 次 在 主 从 群 中 个 具体 问题 的 时 候 都 要 在 解 的 精度 和 计算 效率 之 间 做 折 中 ， 
分 别 执行 全 局 探索 能 力 强 的 CLS 和 开发 能 力 强 的 改进 策略 。 如 果 信 息 交换 的 过 于 频繁 ， 算 法 运行 时 间 会 大 大 增加 ， 效 率 
CDL， 最 后 针对 初始 化 随机 性 过 高 问题 ， 采 用 拉丁 超 立 方 体 低下， 如果 交 换 的 间隔 过 于 长 ， 又 会 使 解 的 质量 大 打折 扣 ， 
采样 对 算法 重建 输入 分 布 ， 在 多 群 和 适应 子 群 职能 更 新 规则 ”在 SEM 中 迁移 周期 的 取 值 设置 为 7。 
双 策略 的 加 持 下， 算法 的 搜索 效率 会 大 大 提升 ， 同 时 从 群 更 (2) 要 交流 的 对 象 。 如 果 在 太 多 群集 之 间 交 流 ， 会 影响 算 
新 规则 的 提出 不 仅 可 以 解决 “维度 退化 ”问题 , 还 避免 了 “过 ”法 的 运行 效率 ， 但 车 交流 的 太 少 ， 又 会 有 很 大 的 陷入 “局 部 
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2.1 


E 细 介绍 CDL-MCPSO 算法 。 
主 从 范式 集群 结构 


为 了 提高 算法 其 


式 的 集群 结构 ， 
( Master ) 和 四 


将 整个 种 群 


探 和 开发 之 间 的 平衡 ， 采 用 
I 分 为 两 类 群体 : 


个 从 群 ( Slave ), EEF 


学 习 策略 在 


略 -综合 维度 


TF YRKE CAR 


到 的 较 优 解 附近 


的 开 


发 ， 主 从 范式 集群 结构 示意 图 如 图 


本 地 最 优 


1 所 示 。 


本 地 最 优 


全 局 最 优 


从 群 2 


全 局 最 优 
从 群 4 


2.2 主 从 群 更 新 策略 


为 使 从 群 能 够 在 局 部 最 优 
改进 后 的 DLS(CDL E ARE 3 


图 1 主 从 范式 集群 结构 
Fig. 1 Master-slave paradigm 


献 [19] 中 被 证 明 


附近 


È 
了 很 强 的 开发 能 力 , 但 是 


发 ”问题 ， 会 极 易 导 致 粒子 陷入 局 部 最 优 ， 


进行 高 精度 的 


通过 


阮 则 ,未 被 改进 的 DLS 在 广 
它 存在 的 “过 度 开 


最 优 ” 的 风险 ， 影 响 解 的 质量 。 同 时 由 于 从 群 主要 负责 在 局 
部 最 优 解 附近 5 进行 高 糖度 的 开发 ， 执行 的 均 为 DLS， 如 果 陷 
基于 主 从 范 。 入 局 部 最 优 ， 利 用 具有 勘探 能 力 的 主 群 搜 索 到 的 信息 即 可 帮助 
一 个 主 群 逃离 ， 所 以 从 税 之 间 无 须 进 行 交流 ， 如 果 交流 反而 会 增加 算 


执行 全 局 搜索 能 力 强 的 综合 


搜索 空间 内 进行 大 范围 的 探索 ， 从 群 执行 改进 策 


法 运行 时 间 ， 所 以 SEM 只 在 主 从 群 之 间 进 行 信息 的 交流 。 
(3) 需 要 交换 的 信息 。 群 集 之 间 通 常 交 换 的 信息 可 以 为 全 


p t 


进行 高 精度 


cluster structure 


PK 


FR, uf 


对 其 改进 : 


去 掉 DLS 速度 更 新 公式 中 对 全 局 最 优 的 学 习 项 , 防止 粒子 个 
几 代 学 习 产 生 的 样本 X? 和 第 三 项 无 限 接近 


kit xr 


经 过 


而 导致 快速 向 全 局 


最 优 Xe%e 靠拢 ， 加 入 CLS 


局 最 优 、 局 部 最 优 以 及 粒子 邻居 最 优 ， 在 CDL-MCPSO 中 ， 
主 从 群 只 在 各 自 迭 代 时 涉及 对 其 他 粒子 的 学 习 ， 协 作 过 程 只 
是 用 来 改善 从 群 可 能 出 现 的 “局 部 最 优 ” 和 主 群 由 于 多 样 性 
强 收敛 速度 不 够 快 的 情形 ， 因 此 无 须 交 换 邻 居 最 优 信息 ， 只 
需 交 换 主 群 全 局 最 优 和 各 从 群 本 地 最 优 信息 。 

(4) 整 合 策略 。 处 理 交 流 后 的 信息 ， 用 来 指导 粒子 的 后 续 
搜索 。 整 个 解 交 换 机 制 的 详细 过 程 总 结 如 下 , 主 群 更 新 策略 : 
在 主 从 群 独立 执行 各 自 的 算法 7 代 之 后 ， 所 有 从 群发 送 自己 
的 本 地 最 优 ibesi(i=1,2,3,4) A ER, 主 群 根据 (1) 在 所 有 从 群 本 
地 最 优 Ibest: G —1,2,3.4) 中 最 优 的 那个 Xis ( sbest )(2) 主 群 自己 
的 本 地 最 优 X%s (3) 主 群 粒 子 自身 个 体 最 优 X%s 对 粒子 进行 
更 新 ， 这 种 更 新 策略 是 基于 主 群 和 从 群 之 间 的 竞争 来 调整 粒 
子 运 行 轨迹 。 速 度 更 新 公式 如 下 : 


Vi G-D- oV (D+ on Gn, (D— X50) 


公式 第 二 项 粒 


Ges (PY 0) XI) (7) 
(1—9)csn, CP; () - X5) 
0, Gbest" < Gbests 
o —10.5,Gbest" = Gbest? (8) 
LGbest" > Gbest? 
其 中 : e 是 协作 因子 ， 位 置 更 新 仍然 采用 式 (2)。 从 群 更 新 策 


各 : 将 所 有 从 群 全 局 最 优 Poss 与 主 群 全 局 最 优 Pies 进行 比较 ， 
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p 


行 了 信息 


Tab. 1 


随机 选择 


RI 解 交 换 机 制 


展示 了 解 交 换 机 制 以 及 多 


Solution exchange mechanism 


上 的 一 个 粒子 ， 
< 享 的 同时 也 不 会 对 算法 的 效率 
E 很 大 影响 ， 所 提出 的 新 的 解 交 换 机 制 不 仅 平衡 了 算法 的 
勘探 和 开发 ， 还 加 快 了 算法 的 收敛 速度 ， 同 时 解 的 质量 也 得 
了 很 大 的 提升 。 表 1 和 表 2 分 别 
f 协 作 相 关 变 量 。 


序 


r1 
号 


具体 操作 


1 
2 
3 
4 
5 
6 


在 每 个 


将 和 ss s. 


所 有 从 群 5S 发 送 自己 的 本 地 最 优 lbesi;(i=1 
计算 Gbests =MIN(f (lbest,(i —1,2,3,4))) 


Ps 按照 式 (7)(2) 对 主 群 更 新 


2,3,4) 给 主 群 M 


计算 Gbest = MIN( Des Pes) 
从 群 中 随机 选 出 一 个 粒子 X 212,34 


今 X,(i 21,2,3,4) = Gbest 


表 2 多 群 协作 相关 变量 
Tab.2 Multi-group collaboration related variables 
变量 变量 含义 
MIS 主 / 从 群 
[m 主 群 全 局 最 优 
Dites 主 群 全 局 最 优 
X Bes 主 群 粒 子 i 个 体 最 优 
Ibest;(i 21,2,3,4) 从 群 本 地 最 优 
Gbest™ 主 群 全 局 最 优 适 应 度 值 
Gbest* 从 群 全 局 最 优 适应 度 值 
Gbest 群体 全 局 最 优 


间隔 划分 ， 
代表 每 个 区 间 的 
较 少 样本 的 前 提 下 
了 在 2 维 搜索 空 
立方 体 采样 产 9 
有 效 避 人 免 算法 


2.4 ”基于 拉丁 超 立方 采样 的 初始 化 策略 


拉丁 超 立 方 


体 抽样 最 早 由 MeKayP^ T. 1 
普通 抽样 方法 蒙特 卡 罗 相 比 ， 它 通过 对 输入 搜索 空间 进行 等 


在 每 个 区 间 中 随机 地 抽取 样本 ， 抽 样 被 强制 地 


IH. 


具有 均匀 分 布 的 特点 ， 


能 够 保 i 
备 比 普 通 样本 更 好 的 高 效 性 ， 


979 年 提出 ， 和 


E 在 拥有 


al 
图 2 展示 


EE mM 


Stepl : 确定 抽 


的 样本 能 更 好 地 均匀 分 布吉 
现 早 熟 收敛 的 现象 。 


间 中 拉丁 超 立 方 体 采样 示意 


FEUX H ; 


图 。 通 过 拉丁 超 
3E 


整个 空间 ， 可 以 
URR FES IRU TF: 


Step2 ; 将 每 维 变量 成 的 定义 域 区 间 Do Xa] 划分 成 五 个 


图 2 和 图 3 分 别 为 拉丁 超 立方 体 采 村 


相等 的 小 区 间 ， 这 样 就 一 共有 no 个 小 超 立方 体 产生 ; 
Step3 ; 产生 一 个 Hxn 的 矩阵 4 ， 
不 2… 好 的 一 个 随机 全 排列 ; 
Steph; A 的 每 行 就 只 有 
小 超 立 方 体内 产生 一 个 样本 ， 这 样 就 共有 五 
选 出 的 样本 互 不 相同 。 


4 的 每 列 均 是 数列 


个 小 超 立方 体 被 选中 ， 在 每 个 


个 样本 被 选 出 ， 
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挑选 出 较 优 的 那个 蔡 换 在 所 有 从 群 中 
这 样 从 群 在 间接 进 
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初始 化 参数 ps niK,m， 基 于 拉丁 超 立方 体 采 样 
产生 初始 种 群 


计算 粒子 适应 度 值 ， 并 技 昭 适 应 度 值 大 小 排序 划分 种 
BÉ (一 个 主 群 ， 四 个 从 群 ) 


初始 化 每 个 粒子 的 速度 V 和 位 置 X 


主 群 (Masten) 执 行 综合 学 习 策略 (CLS) 
| JB Slave rr es ERES ACD) 


5 B E / iag / 
| a / ” 优 位 置 gbest / 
是 mm — 

L————— — sEHRGEENMERIE 


图 3 CDL-MCPSO 算法 流程 
Fig.3 The algorithm flow of CDL-MCPSO 


3  CDL-MCPSO 算法 合理 性 分 析 


CLS 中 去 掉 对 全 局 最 优 项 学 习 操 作 的 合理 性 

由 于 CDL 是 基于 DLS 速度 更 新 公式 开发 ， 因 此 首先 需 
要 对 DLS 速度 更 新 所 采用 策略 存在 的 问题 进行 详细 分 析 , 再 
阐述 为 何 要 进行 CDL 改进 的 原因 ,具体 分 析 如 下 : 和 大 多 数 
PSO 变 体 速度 更 新 公式 一 样 ，DLS 速度 更 新 公式 也 包含 三 个 
部 分 : 惯性 部 分 oV. O 、 个 体 最 优 向 全 局 最 优 学 习 构造 学 习 样 
本 ， 向 样本 学 习 部 分 on; (50 -X,00 、 向 全 局 最 优 学 习 部 分 
cn, XP 0 7 X,,00 ; 本 文 一 般 需要 在 公式 后 两 项 分 别 赋予 种 
群 探索 和 开发 能 力 ，DLS 更 新 公式 最 后 一 项 对 全 局 最 优 
XPO 的 学 习 毫 无 疑问 是 赋予 了 种 群 开 发 能 力 ， 那 么 在 种 群 
已 经 具备 开发 能 力 的 前 提 下 ， 公 式 第 二 项 个 体 最 优 XO" OT 
学 习 对 象 仍然 选择 了 会 进一步 增强 种 群 开 发 能 力 的 全 局 
优 ， 即 使 个 体 最 优 向 全 局 最 优 学 习 后 得 到 的 新 位 置 不 是 
局 最 优 ， 但 是 从 其 具体 的 学 习 过 程 来 看 ， 它 是 逐 维度 进 
学 习 ， 每 一 次 的 学 习 并 不 一 定 只 学 习 一 个 维度 ， 当 得 到 
的 新 位 置 的 适应 度 值 优 于 个 体 最 优 适应 度 值 时 学 习 过 程 
才 会 停止 ， 这 种 学 习 过 程 即 使 放 在 高 维 问题 中 ， 也 会 在 若 
干 代 的 学 习 之 后 使 得 公式 第 二 项 最 终 得 到 的 学 习 样 本 
X1) 变 成 全 局 最 优 , 而 不 是 介 于 个 体 最 优 和 全 局 最 优 之 间 
的 一 个 值 ， 此 时 对 于 后 续 粒 子 的 更 新 ， 公 式 后 两 项 将 会 完 
全 一 样 ， 赋 予 公 式 最 后 一 项 的 开发 能 力 由 于 后 两 项 的 接近 
其 至 完全 一 致 而 变 成 两 倍 ， 若 此 时 粒子 找到 的 全 局 最 优 是 
且 和 迭代 过 程 中 没有 向 其 他 粒子 进行 学 习 ， 
无 疑问 种 群 将 面临 陷入 局 部 最 优 并 且 无 法 逃逸 的 风险 ， 
要 说 明 的 是 ， 公 式 第 二 项 的 学 习 样 本 在 变 成 全 局 最 优 时 
历 的 迭代 次 数 和 本 文 设 置 的 最 大 友 代 次 数 相 比 ， 只 会 占 
小 部 分 比例 。 
针对 以 上 分 析 存 在 的 不 足 ， 本 文 提 出 了 CDL 改进 策略 ， 


3.1 


Zu og AD 


N 


在 二 维 空间 的 示意 


体 原 因 : 首先 为 前 弱 粒 子 在 之 后 大 部 分 迭代 时 期 里 拥有 的 


以 及 所 提 算 法 的 流程 示意 图 。 
图 2 拉丁 超 立 方 体 采样 (2 维 ) 


Fig.2 Latin hypercube sampling(2d) 


两 倍 开 发 能 力 , 需要 去 掉 其 中 一 项 (对 全 局 最 优 的 学 习 ), 同时 
需要 为 种 群 引入 增强 探索 能 力 的 策略 ， 而 针对 多 模 态 问题 6 
发 的 CLS 策略 ， 被 证 明 在 探索 方面 效果 显著 ， 它 能 够 在 含有 
多 
所 


个 峰值 的 优化 问题 中 找到 最 优 解 ， 需 要 说 明 的 是 ， 本 文 之 
以 去 掉 的 是 对 全 局 最 优 的 学 习 项 ， 而 不 是 构造 的 学 习 样 本 
项 ， 是 因为 在 初始 迭代 过 程 中 ， 为 了 不 遗漏 任何 可 能 包含 最 
望 区 域 ， 本 文 希望 算法 在 整体 ] 
发 能 力 ， 而 留 下 来 的 学 习 样 本 项 在 初始 阶段 ， 学 习 样 本 的 值 
还 不 是 全 局 最 优 ， 因 此 对 粒子 的 牵引 力 不 算 很 明显 ， 
可 以 在 搜索 空间 进行 最 大 范围 地 探索 ， 在 后 期 学 习 样 
本 变 为 全 局 最 优 之 后 ， 完 全 可 以 取代 去 掉 的 全 局 最 优 项 来 发 
挥 其 开发 能 力 ， 能 够 使 算法 在 不 同和 迭代 时 期 最 大 程度 地 发 挥 
勘探 或 开发 能 力 ， 综 上 所 述 CDL 策略 中 去 掉 对 全 局 最 优 项 
的 学 习 是 合理 的 。 
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张 其 文 ， 等 : 


3.2 CDL-MCPSO 克服 过 度 开 发 问题 的 能 力 及 理论 性 依据 


CDL-MCPSO 具备 克服 过 度 ] 
主要 有 三 个 : 1、 基于 主 从 范式 的 多 群 策 略 一 能 够 弥补 单个 种 


f 发 问题 能 力 的 理论 性 依据 


B 


f 出 现 过 度 开 发 情 


也 种 群 信息 指导 的 缺点 ，2、 


从 群 的 执行 策略 一 在 


开发 能 力 的 前 提 下 加 入 对 其 他 


粒子 的 学 习 ， 不易 出 现 过 度 开发 ; 3、 解 交换 机 制 (SEM) 的 加 


持 一 即使 开发 过 度 陷入 
作 来 帮助 实现 逃逸 。 


了 最 优 ， 也 可 以 通过 信息 的 交流 协 


蓄 于 主 从 范式 的 多 子 群 相 比 六 


能 够 减 小 粒子 出 现 过 度 开 发 陷入 
以 利用 其 他 种 群 搜索 到 的 信息 来 进行 指导 ， 

EF， 让 主 群 执行 在 多 模 态 
IK (CLS) AIR T E FE HE 


有 更 好 的 多 样 性 


局 部 最 优 的 可 能 性 ， 因 为 可 


其 次 是 在 每 个 子 


群 运行 算法 的 选择 


显 车 的 综合 学 习 策 


行 大 范围 探索 的 能 


问题 中 优化 效果 


力 ， 显 而 易 见 ， 主 群 不 易 出 现 过 度 ] 


BS， 对 于 每 个 从 


群 本 文 执行 优化 后 上 


Im 


赋予 算法 开发 能 力 , 还 


1 


rag 


N 


加 入 了 对 所 有 其 


cr 


ng 


vA 


过 


过 度 开发 的 情形 , 即使 
也 可 以 通过 信息 的 交换 来 增加 种 群 


的 解 交 换 机 制 (SEM)， 
， 而 种 群 多 样 性 世 


LL 
nr 
K 


增加 换 句 话 也 可 以 说 成 是 对 可 


或 一 种 改善 。 


3.8 CDL-MCPSO 不 同和 迭代 时 刻 寻 优 能 
和 开发 能 力 的 平衡 性 
算法 迭代 前 期 ， 将 种 群 划分 为 多 个 子 群 


方 体 采 样 方法 对 种 群 进行 初始 化 ， 保 证 


不 会 使 算法 出 现 


cr 


的 过 度 开 发 情形 的 避免 


的 变化 以 及 勘探 


采用 拉丁 超 立 


种 群 粒 子 在 搜索 空 


间 中 是 均匀 分 布 的 ， 多 样 性 


很 好 的 勘探 能 力 ， 


迭代 过 程 中 主 杂 


FAUT R 


体 最 优 通 过 向 所 有 


也 粒子 的 个 体 最 人 


了 所 有 其 他 粒子 到 


举 本 由 于 是 粒子 的 个 
进行 概率 性 学 习 ， 综 
前 为 止 找到 的 优异 信息 ， 从 群 的 速度 


更 新 公式 中 虽然 力 
法 在 前 期 迁 代 过 和 


比 x? 以 及 主 群 着 
不 明显 ， 能 够 保 记 
算法 迭代 中 期 ， 出 


的 学 习 项 ， 但 是 在 
从 群 学 习 样本 X? 对 粒子 的 指引 作用 相 
期 迭代 过 程 中 发 挥 


> 


H EEG AR HH Ode 


充分 地 勘探 。 
f 中 每 个 粒子 找到 的 个 体 最 优 值 
， 这 是 由 于 迭代 前 期 种 群 


优秀 的 勘探 能 力 所 发 挥 的 作用 而 产生 的 结果 ， 此 时 不 论 是 对 
于 主 群 还 是 从 群 学 习 检 于 被 学 习 对 象 的 解 变 优 ， 
使 得 学 习 样本 的 值 也 随 之 变 好 ， 整 个 种 群 的 更 新 都 在 朝 最 终 
的 最 优 解 慢 慢 靠近 ， 开发 能 力 并 存 ， 优 化 了 解 的 
质量 但 不 会 使 粒子 过 勘探 能 力 能 够 继 
续 探 索 可 能 的 较 优 解 。 
SHE TUS BI. E 9a S BE ds 
从 群 中 学 习 样本 X 的 值 越 来 越 接近 直至 等 于 全 局 最 优 , 而 另 
外 的 学 习 样 本 xr 即 使 是 向 其 他 粒子 的 个 体 最 优 进行 学 习 所 
得 ， 但 由 于 迭代 过 程 已 经 至 的 后 期 ， 此 时 每 
粒子 找到 的 个 体 最 优 解 的 太 ] 题 的 最 终 解 很 近 了 
时 Xr 的 值 浮动 逐渐 减 小 了 进行 大 范围 的 探索 
要 围绕 在 当前 找到 的 最 优化 ， 因 为 前 期 展 好 
勘探 能 力 和 中 期 适应 的 勘探 能 力 使 得 粒子 对 问题 空间 已 经 


高 度 地 探索 ,所 以 在 后 全 局 最 优 的 过 程 中 ， 


EEREZEEG 
EE 
ex 


RE 


不 存在 靠近 的 是 
的 是 开发 能 


综 上 ， 改 进 粒子 群 算法 的 寻 


代 时 刻 的 不 同 需求 分 另 
和 开发 能 力 ， 实 现 了 两 者 之 间 昌 
3.4 CDL-MCPSO 算法 的 收敛 性 分 析 


亿 能 力 能 够 根据 算法 不 同 迭 
I 发 挥 勘 探 或 开发 能 


此 在 后 期 整个 种 群 主 


或 同时 发 挥 勘 探 


Solis 等 人 P 的 研究 结果 证 明了 随机 优化 算法 以 概率 1 


收敛 于 全 局 最 优 解 的 条 伯 
假设 1 若 f(DG,O))<f09， 


T 


机 


1 和 1 


设 2， 就 可 以 保 订 
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f(D(x,6)) < f (0) 
标 函数 ，D 为 产生 问题 解 的 函数 ，x 为 5 
中 搜索 到 的 一 个 点 ，“ 为 从 概率 空间 GB. 产生 的 随 
5 为 问题 的 约束 空间 。 


S 的 任意 Borel TÈ A, 若 其 测度 v0 >0， 则 有 : 


I[a-4u»-0 
p=] 


h，vC >o0 为 子 集 4 D n Hf 
的 概率 。 


EHE, uA 为 由 测度 产生 4 


假设 目标 函数 f 为 可 测 函 数 ， 问 题 的 约束 空间 
5 为 可 测 子 集 , {3% 为 随机 算法 产生 的 解 序列 , 则 当 满 足 假设 


段 设 2 时 ， 有 : 


limP[x € R;,]-1 


R 为 全 局 最 优点 的 集合 。 
此 ， 对 任 一 随机 优化 算法 ， 只 要 能 够 满足 假设 1 与 假 


DG, X) -| 


FE 它 能 以 概率 1 收敛 于 全 局 最 优 解 。 

定理 2 CDL-MCPSO 算法 能 够 以 概率 1 收敛 于 全 局 最 
优 解 

证 明 CDL-MCPSO 算法 满足 假设 1。 

设 CDL-MCPSO 算法 的 迭代 


函数 DD 可 以 表示 为 


Pero F0 S f (xi) 
Xiro f Xa) Sf (Poi) 


E 明 CDL-MCPSO } 


Riz 1. 


CDL-MCPSO 算法 满足 假设 2。 


设 规模 为 N 的 粒子 群 的 样本 空间 的 并 必须 包含 $5， 即 
P, Mu KPRF i 的 样本 空间 的 支撑 集 。 对 于 


满足 %'=P=Ps 的 粒子，Mis=5 ， 其 他 粒子 i 则 满足 : 


对 标准 PSO 9 
人 入 位 置 更 新 公式 后 ， 代 入 Mu 可 得 : 


ERESZ] <a, 0<, Si s 


Ri H A QR s (9-—C0n, 将 速度 更 新 公式 


Mi, = Xij t OX jaa = Xi ja-2) T 


P (Pija —XÀ ja) 9s (m — Xi ja ) 


I ERA 


显然 ，Mi, 表 示 一 个 由 9 p ME 


E 形 , 其 中 一 个 端点 为 9 =p=0, 男 一 个 为 9=0 ,=c。 


当 max(c, [Pia — Xi aal € [Pye = Xi ja- p «0.5x diam;(S) 成 立时 , 显 


AUR vM u(19 «v6 。 其 中 ，diam,(5) RIR S 在 第 J 维 分 量 的 长 


度 。 
Aa 


HT 
生 新 的 粒子 。 
以 上 改进 策略 产生 


C ， 


NH, 82S. 


s, 


于 x, (ap; * €) p,)/ (c +c) š 
的 增加 ， VM) 不 断 减少 ， 
(Que, 表示 存在 整数 t, 使 得 当 t»t, 存在 集合 ACS, 使 


局 最 优 解 。 


CDL-MCPSO 4i f fi 


于 改进 策 


RI 


值 的 评 


十 以 及 速度 夭 


Pg mM 0, 从 而 ， 随 进 代 次 


H 并 Quo 也 在 减少 。 从 而 


得 jw-*。， 这 表明 标准 PSO 算法 不 满足 假设 2， 不 能 以 概率 


1 收敛 了 


FE 标准 PSO 算法 的 基础 上 ， 运 
# 、 不 同学 习 策 略 以 及 种 群 之 间 的 协作 交换 机 制 ， 
丸 此 ， 对 正常 进化 的 粒子 ， 设 其 支撑 集 的 并 集 为 
的 粒子 ， 设 其 支撑 集 的 并 集 为 8 。 
各 的 指导 性 ， 必 然 存在 整数 点， 使 得 当 z>a 
因此 对 于 CDL-MCPSO 算法 , 存在 整数 ,使 得 
f, eUZ35., 
定义 5 的 任意 Borel FH A7 M. , WA A 90 , sa-a, 
证 -ww-o， 所 以 CDL-MCPSO 算法 满足 假设 2. 


CDL-MCPSO 算法 是 全 局 收敛 算法 。 


Wing CDL-MCPSO 算法 生成 的 解 序列 ， 
于 CDL-MCPSO 算法 满足 假设 1 和 假设 2， 通过 定理 可 知 : 
limPp, sR&]=1I 成 立 ， 故 CDL-MCPSO 算法 能 够 以 概率 1 收敛 
于 全 局 最 优 解 。 

算法 的 收敛 速度 和 解 的 精度 在 实验 部 分 进行 了 详细 分 析 ， 见 后 文 。 
3.5 算法 的 时 间 复 杂 度 分 析 

传统 PSO 算法 主要 包括 以 下 几 个 部 分 : 初始 化 、 适 应 度 
I 位 置 更 新 ， 他 们 的 复杂 度 分 别 为 


录用 定稿 张 其 
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O(QxNxD). O(NxD)AIl O(2XNxD)(N 和 D 分 别 表示 种 群 规 模 
和 维度 。 因 此 ，PSO 的 时 间 复 杂 度 为 ONxD)。 与 传统 PSO 


相 比 ，CDL-MCPSO 算法 需要 构建 学 习 样本 、 主 从 群 信息 交 


流 时 评估 所 有 从 群 本 地 最 优 。 由 于 主 从 群 只 


在 连续 和 迭代 次 


式 第 二 项 得 到 的 样本 已 经 接近 全 局 最 优 ， 但 主 群 中 其 他 粒子 


个 体 最 优 的 差异 仍然 保持 学 习 后 的 样本 
DLS 的 一 部 分 开发 和 


就 会 进行 解 


言 息 的 交换 ,粒子 只 是 慢 慢 开始 出 现 
但 仍 具备 较 好 的 种 群 多 样 性 。 


数 达 到 刷新 闻 阶 时 才 进 行 信息 交流 ， 并 且 评 估 的 本 地 最 优 分 
别 来 自 四 个 从 群 ， 每 次 交流 只 评估 四 组 数据 ， 因 此 多 群 协作 
过 程 的 时 间 复 杂 度 为 0(3)， 并 且 只 有 当 粒 子 的 个 人 最 佳 位 置 
更 新 或 个 人 最 佳 位 置 在 儿 代 内 没有 更 新 时 ， 才 需要 重建 学 习 
样本 。 构 建 学 习 样本 的 最 坏 情 况 时 间 复 杂 度 是 OND), 判断 


迭代 是 否 达 到 刷新 间隙 以 及 是 否 停 
Xl, CDL-MCPSO 的 最 坏 1 


根 
基本 PSO 算法 处 于 同一 水 平 。 


止 的 时 间 复 杂 度 为 O(2), 
青 况 时 间 复 杂 度 也 是 OND), 包 
括 初 始 化 O(2xNxD), 评 估 O(NxD+5), 更 新 O(2xNxD+NxD)。 
居 以 上 分 析 ， 可 以 得 知 CDL-MCPSO 算法 的 时 间 复 杂 度 与 


备 较 好 的 多 样 性 ， 
E 力 被 中 和 ， 加 上 种 群 每 隔 一 定 迭 代 次 数 
改 敛 的 趋势 ， 


6 代表 种 群 在 迭代 后 期 的 大 致 分 布 ， 从 群 迭代 
DLS 公式 第 二 项 的 学 习 样 本 等 于 全 局 最 优 ， 后 两 项 均 发 
挥 开 发 作用 ， 种 群 多 样 性 减弱 ， 大 部 分 粒子 聚集 在 全 局 
最 优 解 附近 , 但 CLS 的 存在 , 仍然 会 使 个 别 粒子 散落 在 
距离 全 局 最 优 解 不 太 远 的 位 置 。 
4.2 测试 函数 
为 验证 CDL-MCPSO 性 能 , 采用 Benchmark 测试 函数 集 
中 的 10 个 函数 进行 求解 验证 , 其 中 包括 两 个 单 峰 函数 f — f ~ 
五 个 多 峰 函 数 请 - 方 、 三 个 组 合 函数 -fo， 测 试 函数 的 县 体 
Hd ülxe3 所 示 。 
表 3 实验 使 用 的 10 个 测试 函数 
Tab.3 The 10 test functions used in the experiment 
No. Functions FE =F(x’) 
Rotated High Conditioned Elliptic 
Unimodal f . -1300 
. Function 
Functions 
h Rotated Discus Function -1100 
f Rotated Weierstrass Function -600 
] fa Rotated Griewanks Function -500 
Basic 
. Ís Rotated Katsuura Function 200 
Multimodal . D. . 
) fs Lunacek Bi Rastrigin Function 300 
Functions 
Rotated Lunacek Bi Rastrigin 
f . 400 
Function 
Composition Function 2 
fs 800 
(n-3,Unrotated) 
Composition Composition Function 3 900 
Functions ° (n-3,Rotated) 
Composition Function 4 
fio 1000 


(n=3,Rotated) 


4.3 ”对比 算法 


为 了 验证 CDL-MCPSO 算法 在 求解 复杂 问题 时 的 性 能 ， 


将 其 


与 统一 粒子 群 算法 (unified particle swarm optimization ， 


UPSO)U03、 完 全 知情 粒子 群 算法 (fully informed particle swarm, 
FIPS)00 维 度 学 习 粒 子 群 算法 (particle swarm optimization 
based on dimensional learning strategy，DLPSO)20、 并 行 综合 


4 ”实验 验证 与 分 析 
4.1 典型 迄 代 时 刻 的 个 体 分 布 图 及 种 群 多 样 性 的 分 析 
给 出 了 种 群 分 别 在 50、500 和 900 代 时 的 个 体 分 布 图 ， 
如 图 4-6 所 示 。 
PES 
2f u "E 
" " 口 二 .= og cH 口 口 
EE o ge 7 : 
E o LI 
Wa = n a 口 
3 口 bb i 
2 o 
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粒子 序号 
图 4 ”种群 在 50 代 时 的 个 体 分 布 图 
Fig.4 Individual distribution of the population at the 50th generation 
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图 5 种 群 在 500 代 时 的 个 体 分 布 
Fig.5 Individual distribution of the population at the 500th generation 
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图 6 种 群 在 900 代 时 的 个 体 分 布 
Fig.6 Individual distribution of the population at the 900th generation 
图 4 代表 种 群 在 迭代 前 期 的 大 致 分 布 ， 主 从 群 独立 执行 
CLS 和 DLS，CLS 学 习 对 象 为 所 有 其 他 粒子 的 个 体 最 优 ， 


DLS 公式 第 二 项 学 习 对 象 虽 为 全 局 最 优 ,但 得 到 的 样本 距离 
还 要 些许 途 代 时 间 ， 种 群 多 样 性 良好 ， 粒 子 基 


成 为 全 局 最 
本 均匀 散落 在 搜索 空间 


图 5 代表 种 群 在 迭代 中 期 的 大 致 分 布 , 从 群 闪 代 DLS A 


学 习 粒 子 群 优化 算法 (parallel comprehensive learning 


particle 


swarm optimization, PCLPSO)07 及 标准 粒子 群 算法 (A 


modified particle swarm optimizer, PSO-W) 
算法 进行 对 比 实验 。 为 了 保证 测试 的 公 习 


IN 


五 种 改进 粒子 群 


性 ， 算 法 的 参数 设 


置 均 相 同 : 种 群 规模 为 40,， 迁 代 次 数 为 1000， 所 有 算法 分 别 
在 D-10. 30. 50 三 个 维度 下 在 每 个 测试 函数 上 独立 运行 30 
次 。 实验 环 境 设 置 如 下 : AMD Ryzen 5350H with Radeon Vega 
Mobile Gfx 2.10 GHz，RAM 16GB，Windows 10 操作 系统 ， 
Matlab R2016b。CDL-MCPSO 算法 及 对 比 算法 的 参数 设置 如 表 4 
所 示 。 


表 4 CDL-MCPSO 算法 及 对 比 算法 的 参数 设置 
Tab.4 CDL-MCPSO algorithm and comparison algorithm parameter setting 


算法 参数 设置 参考 文献 
PSO-W a=0.9~0.4,c1=c2=2 [2] 
UPSO 3:0.729,c71.49445 [13] 
FIPS X:0.729,c=2 [11] 
DLPSO w=0.9~0.4,c1=1.5,c2=0.5~2.5 [21] 
PCLPSO w=0.9~0.4,c=1.49445 [17] 


CDL-MCPSO @w=0.9~0.4,c1=c2=2.05,c3=2.0,c=1.49445 
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4.4 算法 性 能 分 析 函数 、 多 模 态 函数 还 是 组 合 函数 上 均 能 找到 更 好 的 或 相当 的 优 
dé 5-6 表示 在 D=30,50 两 种 维度 情况 下 ， 每 种 算法 在 测试 MWER, ie 在 个 别 函 数 上 稳定 性 稍 显 不 足 ， m Qi 
函数 上 的 平均 值 和 标准 偏差 。 六 种 算法 的 最 佳 结果 均 以 粗 体 显 ”的 主要 原因 是 主 从 群 采用 的 综合 学 习 策 略 (CLS) 和 综合 维度 学 
示 出 来 。 从 表格 数据 可 以 看 出 CDL-MCPSO 算法 除了 在 测试 函 习 策 略 (CDL) 在 算法 迭代 前 期 发 挥 了 优秀 的 勘探 能 力 ,探索 了 可 
数 户 上 的 性 能 略 显 不 足以 外 ， 能 够 在 其 他 绝 大 多 数 不 管 是 单 峰 ”能 存在 较 优 解 的 所 有 区 域 ， 为 最 终 解 质量 的 提高 莫 定 了 基础 。 


表 5 对 比 算法 在 测试 函数 上 的 优化 结果 (D=30) 


Tab.5 The optimization results of the comparison algorithm on the test function(D=30) 


函数 Criterion PSO-W UPSO FIPS DLPSO PCLPSO CDL-MCPSO 

Mean 1.07E+08 1.70E+07 1.67E+08 5.36E+07 3.26E+08 1.60E+07 

fi Std 8.06E+07 1.38E+07 1.86E+07 6.23E+07 8.16E*07 4.55E+06 
f. Mean 1.04E+05 5.45E+04 8.45E+04 6.27E+04 3.95E+04 4.79E+04 
Std 3.01E+04 1.67E--04 2.22E*04 2.04E*04 2.26E*04 4.10E+04 

Mean 2.52E+02 2.47E+02 3.93E+01 3.12E+01 3.06E+01 2.39E+01 

" Std 1.S1E+00 3.25E+02 3.78E+00 1.74E+00 3.35E+00 3.22E+01 
Mean 2.42E+03 4.08E+03 5.78E+03 4.74E+03 4.24E+03 3.76E+03 

5 Std 3.42E+03 2.24E+03 7.14E+02 3.0SE+02 6.18E+02 7.50E+02 
Mean 2.28E+02 2.28E+02 3.09E+00 2.47E+00 2.32E+00 1.33E+00 

Std 3.38E+02 3.39E+02 3.01E-01 4.64E-01 5.96E-01 4.39E-01 
Mean 6.92E*02 3.32E*02 2.35E+02 2.92E+02 1.93E+02 9.32E+01 

f Std 1.26E+02 2.68E+02 2.74E+01 6.07E+01 7.38E+01 1.12E+01 
Mean 6.84E+02 3.83E+02 2.46E+02 2.61E+02 2.39E+02 1.34E+02 

ý Std 1.19E+02 2.34E+02 1.91E+01 5.46E+01 1.05E+02 1.46E+01 
Mean 5.51E+03 3.25E+03 6.74E+03 5.89E+03 3.58E+03 4.29E+02 

á Std 3.49E+03 1.89E+03 2.19E+02 1.56E+03 8.48E+02 3.85E+02 
Mean 8.48E+03 7.77E+03 8.97E+03 7.95E+03 6.61E+03 5.30E--03 

h Std 3.47E+03 2.57E+03 2.52E402 8.83E+02 8.85E+02 2.66E+02 
Mean 4.31E+02 4.21E+02 3.03E+02 2.88E+02 2.94E+02 2.1SE+02 

Std 3.68E+00 2.10E+02 2.04E+02 1.08E+01 7.63E+00 5.01E+00 


表 6 对 比 算法 在 测试 函数 上 的 优化 结果 (D=50) 


Tab.6 The optimization results of the comparison algorithm on the test function(D=50) 


函数 Criterion PSO-W UPSO FIPS DLPSO PCLPSO CDL-MCPSO 
CO Mean 6.63E+08 4.03E+08 3.64E+08 2.12E+08 2.27E+08 2.60E+07 
pum 1 Std 2.68E+08 3.10E+07 5.13E+07 3.74E+07 1.64E+08 1.0SE+07 
l f Mean 9.46E+04 1.63E+05 1.72E+05 9.73E+04 8.79E+04 9.93E+04 
Std 7.02E+04 6.54E+04 2.03E+04 2.62E+04 2.12E+04 3.42E+04 
Mean 2.75E+02 2.67E+02 7.43E+01 5.76E+01 6.11E+01 4.90E+01 
^ Std 2.27E*00 3.12E+02 3.44E+00 3.60E+00 5.09E+00 2.06E+00 
Mean 1.88E+04 2.17E+04 1.92E+04 1.71E+04 1.90E+04 1.52E+04 
^ Std 6.77E+03 4.36E+02 1.59E+03 1.10E+03 1.67E+03 4.94E+03 
Mean 2.28E+02 2.28E+02 4.15E+00 3.71E+00 3.65E+00 2.21E+00 
i Std 3.37E-01 3.38E+02 3.38E+02 5.89E-01 4.75E-01 1.06E+00 
Mean 1.61E+03 4.54E+02 5.32E+02 6.87E+02 6.81E+02 2.32E+02 
fs Std 7.17E+02 1.92E+02 4.35E+01 1.05E+02 2.98E+02 2.35E+01 
Mean 1.62E+03 5.47E+02 5.30E+02 7.20E+02 6.64E+02 3.76E+02 
n Std 7.21E*02 1.39E+02 2.69E+01 1.11E+02 2.35E+02 9.88E+01 
Mean 1.00E+04 6.80E+03 1.34E+04 1.17E+04 9.00E+03 5.70E+03 
P Std 4.94E+02 4.47E+03 2.08E+03 1.74E+03 1.32E+03 4.75E+02 
Mean 7.68E+03 1.87E+03 1.51E+03 3.23E+03 2.97E+03 1.40E+03 
^ Std 1.93E+03 1.43E+03 1.32E+03 1.75E+03 5.26E+03 1.36E--03 
Mean 4.90E+02 4.74E+02 3.92E+02 3.72E+02 3.77E+02 3.39E+02 
fi Std 2.79E+00 1.77E+02 1.36E+01 1.31E+01 1.57E+01 1.68E+02 


图 7-12 分 别 展示 了 CDL-MCPSO 算法 在 A (HIE Aa (E fi, mi PCLPSO 不 仅 没有 找到 全 局 最 优 解 , 而 且 从 50 代 左 右 就 
模 态 ) 和 (组合) 三 个 函数 上 30、50 维 下 的 收敛 曲线 图 , 针对 函 始 收敛 ,收敛 速度 过 于 快 ， 图 10 和 图 7 的 情况 大 致 类 似 。 针 
数 fi 的 收敛 曲线 图 , 从 图 7 可 以 观察 到 改进 算法 明显 优 于 PSO- ”对 函数 的 收敛 曲线 图 ， 尤其 是 图 11，CDL-MCPSO 性 能 明显 
W, UPSO, FIPS, 5j PCLPSO 相 比 ， 改 进 算法 找到 了 全 局 最 优优 于 其 他 对 比 算 法 ， 从 fs 的 收敛 曲线 图 可 以 直观 地 看 出 CDL- 
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MCPSO 不 仅 在 找到 的 解 上 明显 优 于 其 他 对 比 算 法 , 而 且 收 敛 速 
度 也 快 于 其 他 算法 , 在 迭代 后 期 没有 出 现 “ 振 荡 ” 现 象 , 解决 了 
“维度 退化 ”问题 ， 虽 然 CDL-MCPSO 算法 也 存在 不 足 ， 算 法 


稳定 性 在 个 别 测试 函数 上 相 比 其 他 算法 略 显 不 足 , 但 总 的 来 说 ， 


CDL-MCPSO 算法 与 其 他 算法 相 比 , 在 优化 结果 上 都 大 有 提升 。 


O 100 200 300 400 500 600 700 800 900 1000 
点 代 次 数 


7 KZ 的 收敛 曲线 (D=30) 
Fig.7 Convergence curve of fı (D-30) 


对 


LPSO 
——— CDL-MCPSO 


500 600 700 800 900 1000 
代 次 数 


图 8 函数 的 收敛 曲线 (D=30) 
Fig.8 Convergence curve of f, (D-30) 
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Edo 函数 广 的 收敛 曲线 (D=30) 
Fig.9 Convergence curve of f; (D=30) 
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图 10 函数 fi 的 收敛 曲线 (D=50) 
Fig. 10 Convergence curve of f, (D=50) 


适应 度 依 
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图 11 X 的 收敛 曲线 (D=50) 
Fig. 11 Convergence curve of f; (D-50) 
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函数 fs 的 收敛 曲线 (D=50) 


12 
Fig. 12 Convergence curve of f; (D=50) 
为 了 进一步 验证 CDL-MCPSO 算法 的 性 能 , 借鉴 数据 统 


计 与 分 析 方 法 ， 采 用 显著 性 水 平 为 0.05 的 Wilcoxon 秩 检 验 
方法 来 判断 算法 性 能 。 其 中 “+”“ 一 ”“ 守 ”分 别 表示 CDL- 


MCPSO 算法 的 结果 优 于 、 


果 。 


劣 于 、 相 当 于 对 应 算法 的 测试 结 
从 表 7 的 Wilcoxon 测试 结果 来 看 ， 在 <=0.05 时 CDL- 


MCPSO 算法 在 测试 函数 上 相 较 于 对 比 算法 均 有 一 定 优势 。 


表 7 通过 Wilcoxon 的 测试 得 到 结果 
Tab.7 Results obtained through Wilcoxon's test 


D CDL-MCPSO VS P-value t ~  —  Qg-0.05 
PSO-W 0.000851 37 1 0 Yes 
UPSO 0.000112 36 1 0 Yes 
FIPS 0.000401 36 1 0 Yes 
xi DLPSO 0.000102 35 2 0 Yes 
PCLPSO 0.000099 36 1 0 Yes 
PSO-W 0.000134 36 1 0 Yes 
UPSO 0.000026 35 1 0 Yes 
FIPS 0.000147 37 1 0 Yes 
" DLPSO 0.000082 36 2 0 Yes 
PCLPSO 0.000096 36 1 0 Yes 


HR 


根据 以 上 对 比分 析 ，CDL-MCPSO 可 以 保持 更 好 的 多 样 
F 发 原始 空间 来 获得 高 质量 的 解 ， 并 且 有 具备 最 佳 的 达到 最 


优 的 能 力 使 算法 快速 收敛 。 因 此 ，CDL-MCPSO 是 一 种 能 够 


有 效 提 升 PSO 性 能 的 可 行 方法 。 


5 


结束 语 


为 提高 算法 运行 效率 、 平 衡 算 法 的 勘探 和 开发 能 力 ， 本 


文 基于 主 从 范式 集群 结构 提出 了 综合 维度 学 习 的 多 和 群 协作 粒 
子 群 优化 算法 (CDL-MCPSO), 主 从 群 独立 执行 具有 不 同 职能 


的 学 习 策 略 ， 其 中 从 群 学 习 策 略 是 针对 DLS 存在 的 “过 度 开 


发 ”问题 而 开发 ， 同 时 为 提升 解 的 质量 ， 还 提出 了 一 种 新 的 
解 交 换 机 制 , 用 于 群体 之 间 的 信息 共享 , 实验 结果 证 明 CDL- 


MCPSO 算法 不 仅 能 够 在 绝 大 多 数 单 峰 、 多 峰 和 组 合 函数 上 


找到 更 好 的 或 相同 的 优化 结果 ， 而 且 在 收敛 速度 和 解 的 质量 
上 均 有 明显 优势 ， 不 仅 解决 了 “振荡 ”和 “维度 退化 ”问题 ， 
还 有 效 避 免 了 “过 度 开 发 ” 下 一 阶段 工作 主要 是 提高 算法 在 


部 分 测试 函数 上 的 稳定 性 并 将 其 应 用 于 实际 复杂 问题 中 进行 


有 效 性 的 验证 。 
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